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Dynamic  Plastic  Instabilities  in  Nonlinear  Inelastic  Response  to  Pulse  Loading 
Final  Report 

ARO  Proposal  Number  24362-EG,  Coiitract/Grant  No.  DAAL03-87-K-0038 
(15  March  1987  -  14  September  1991) 

Abstract.  The  report  summarizes  the  main  results  obtained,  with  emphasis  on  the  energy 
approach  introduced  first  for  a  single  degree  of  freedom  beam  model,  and  then  develops  for  one 
of  two  degrees  of  freedom.  The  elastic  strain  energy  function  is  plott^  as  a  surface  over  the 
di:^lacement  coordinate  plane,  whose  topogr^hy  depends  on  the  plastic  strains.  It  provides  a 
guide  to  the  various  types  of  beam  response,  quasiperiodic  and  chaotic,  and  enables  their 
dependence  on  parameters  such  as  impulse  magnitude  and  damping  ratio  to  be  clearly  undmtood. 

1.  Statement  of  problem 

The  subject  of  the  research  project  was  defined  as  the  investigation  of  anomalous  inelastic 
responses  of  certain  structural  elements  (beams  and  plates)  to  short  pulse  transverse  loading.  These 
"anomalous"  behaviors  are  obsaved  in  computed  final  (permanent)  deflections  that  exhibit  extreme 
sensitivities  to  parameters.  In  certain  circumstances,  b^use  of  these  not  even  the  sign  of  the  major 
deflection  can  be  reliably  predicted,  either  for  a  physical  laboratory  specimen  or  "on  paper"  for  a 
mathematically  well  defined  problem. 

These  phenomena  indicate  the  presence  of  multiple  equilibrium  states,  stable  and  unstable,  in  the 
new  geometry  of  the  structure  created  by  plastic  deformations  resulting  tfom  the  pulse  loading. 

Tlie  structures  concerned  are  beams  and  plates  fixed  at  their  boundaries,  so  that  transverse  loading 
which  causes  small  plastic  deformations  effectively  converts  the  beam/plate  to  a  shallow  arch/shell. 
The  dynamic  instabilities  are  thus  related  to  those  involved  in  snap-buckling.  The  fundamental 
difference  is  that  there  the  structure  and  loading  are  known  and  specified  in  advance,  while  in  the 
present  class  of  problems  the  relevant  structure  and  loading  cannot  be  specified  except  by  means  of 
calculated  response  quantities;  plastic  deformations  and  available  energies,  in  particular. 

2.  Summary  of  main  results 

Following  the  accidental  discovery  of  these  anomalous  behaviors  in  1984,  preliminary  research 
received  support  from  an  ARO  "Feasibility  Study"  grant  (Proposal  No.  23225-EG).  Papers  (1-4] 
in  the  list  of  publications  present  results  obtained  prior  to  the  present  grant. 

These  preliminary  results  further  displayed  the  sensitivity  of  the  numerical  solutions  for 
discretized  uniform  beams  to  the  parameters  of  the  problem  and  the  numerical  treatment.  The  IDoF 
model  studies  with  damping  showed  that  the  previous  calculations  for  an  undamped  model  are 
inc^able  of  predicting  the  sign  of  the  final  rest  deflection  of  the  damped  specimen.  The 
"characteristic  diagram"  for  toe  damped  model,  in  which  the  final  deflection  is  plotted  against  toe 
load  pulse  parameter,  showed  the  final  deflection  alternating  abruptly  between  positive  and  negative 
values.  The  discontinuous  transitions  correspond  to  singular  points  in  the  phase  diagram.  TTie 
pattern  of  final  displacements  has  some  regularity,  but  the  dependence  on  damping  straigto  and 
other  parameters  appears  to  be  quite  complicated  [4]. 

The  main  result  of  the  research  supported  by  the  present  grant  is  a  new  approach,  introduced  in 
[5],  to  displaying  and  interpreting  the  response  features,  simpler  and  more  dhwt  than  the  usual 
stability  studies  in  taros  of  phase  space  geometry.  The  discontinuous  alternation  of  the  final  state 
between  positive  and  negative  values,  and  all  other  features  of  the  characteristic  diagram  of  this 
model,  b^me  immediately  understandable  from  energy  diagrams,  constructed  by  plotting  the 
elastic  strain  energy  and  the  total  energy  against  the  displacement,  as  furnished  by  the  numerical 
solution  of  the  equations  of  dynamics  auid  plasticity.  The  approach  is  based  simply  on  the  energy 
balance  equation,  i.e.  the  first  Integra,  of  the  equations  of  motion.  The  sum  of  toe  current  kinetic 
and  elastic  strain  energies  is  equal  to  toe  total  available  energy,  which  can  be  written  as  a  reference 
value  less  the  energies  subsequaitly  dissipated  in  damping  and  in  plastic  work,  if  any.  Applied  to 
the  1  DoF  model,  the  elastic  strain  energy  function  is  a  quartic  polynomial  function  of  toe 
displacement,  with  constant  coefficients  as  soon  as  the  plastic  strains  have  become  constant. 
Typically  it  then  exhibits  a  local  maximum  between  two  local  minima,  i.e.  a  potential  hUl  between 
two  wells.  The  final  rest  state  is  always  at  one  or  the  other  of  the  local  minima. 


An  early  application  of  the  insights  provided  by  these  diagrams  was  to  explain  the  wigin  of 
erroneous  results  in  a  contempoip^  paper  making  use  of  the  same  model.  In  this  paper  the  authors 
considered  loading  specified  by  initial  conditions  (as  in  our  work),  as  well  as  periodic  loading.  In 
the  former  case  they  obtained  a  diagram  in  which  the  coordinates  are  the  initial  velocity  and 
displacement,  and  a  dot  is  entered  if  tlie  final  displacement  is  negative  (opposite  to  the  irr  Lai 
impulse).  Their  diagfM  is  extremely  complex,  with  data  points  suggestive  of  fractal  btnindaries 
between  attracting  Imsins.  The  actual  response  has  no  such  complications;  no  ftactai  geometry  or 
other  signs  of  chaotic  behavior.  It  is  obvious  from  the  energy  plots  tliat  the  numericsai  r  oludon  is 
everywhere  calculable  without  difficulty,  contrary  to  the  authors'  statement  that  aiiempled  numerical 
solutions  are  "meaningless".  The  basic  error  of  the  published  paper  was  explained  [6]  with 
further  development  in  (7].  Numerical  difficulties  are  encountei^  only  if  the  damping  coefficient  is 
taken  excessively  small. 

Publications  [8, 9]  provide  further  examples  of  finite  element  solutions  for  pin-ended  and  fully 
clamped  beams,  illustrating  unexpected  (apparently  chaotic)  features  of  time  histories,  contrasting 
with  tlie  comparative  simplicity  of  the  IDoF  model.  The  energy  approach  was  reviewed  in  these 
papers  for  the  IDoF  model,  and  in  [9]  was  applied  to  cases  of  uniform  beams,  taking  data  from  the 
finite  element  solutions,  and  plotting  the  two  energies  as  functions  of  the  midpoint  displacement. 
These  diagrams  iUustrate  the  expect  main  features  of  damped  elastic-plastic  responses,  but  are 
irregular  and  uninstructive  since  the  strain  energy  is  a  function  not  only  of  the  midpoint 
displacement,  but  of  ail  the  kinematic  variables  of  the  discretized  beam. 

Since  September  1989  the  work  has  concentrated  on  studying  a  two  degree  of  freedom  (2DoF) 
beam  model.  We  again  took  this  to  be  of  Shanley  type,  with  deformable  sandwich  beam  elements 
at  the  two  quarterpoints  and  midpoint.  Assuming  symmetrical  deformations,  the  two  transverse 
displacements  at  ^e  midpoint  and  quarterpoint  specify  the  configuration.  The  phase  space 
(displacements  and  velocities)  now  having  four  dimensions,  a  much  richer  variety  of  behaviors, 
including  chaotic  vibrations,  is  exhibited.  The  use  of  a  2DoF  model  is  advantageous  especially  for 
displaying  the  energy  approach.  The  strain  energy  is  a  function  of  the  two  displacements,  with  the 
four  plastic  strains  appearing  as  psuameters.  When  they  become  constant,  the  function  can  be 
represented  as  a  surface  over  the  displacement  coordinate  plane.  The  important  typical  features  of 
tfiis  surface  created  by  plastic  straining  are  a  local  maximum,  two  saddle  points,  and  two  local 
minima  The  character  of  the  response,  after  the  completion  of  plastic  deformation,  depends  on  the 
position  of  the  total  energy  plane  relative  to  these  five  singular  points.  The  total  energy 
monotonically  decreases,  rapidly  when  plastic  flow  occurs  and  more  gradually  due  to  damping 
losses.  Thus  the  corresponding  plane  steadily  descends,  with  intersection  curves  of  various  shapes 
and  steadily  decreasing  area  The  actual  state  point  must  lie  in  the  interior  of  the  current  intersection 
curve.  If  the  total  energy  plane  is  initially  above  the  top  of  the  potential  hill  (local  maximum),  the 
response  may  be  either  quasiperiodic  or  weakly  chaotic  until  it  reaches  the  vicinity  of  the  maximum 
and  saddle  points.  Then  it  typically  becomes  strongly  chaotic,  with  a  highly  irregular  mixture  of 
large  amplitude  vibrations  and  small  amplitude  positive  and  negative  vibrations.  Eventually  there  is 
an  abrupt  transition  to  a  small  amplitude  vibration,  which  retains  its  sign  and  becomes  smoother 
(quasip^odic)  as  the  amplitude  decreases  to  zero.  This  transition  indicates  that  the  total  energy  has 
f^len  below  the  energy  at  the  lower  saddle  point,  so  that  the  solution  path  is  trapped  in  one  or  the 
other  of  the  two  potential  wells.  The  final  state,  of  course,  is  at  one  of  the  two  local  minima.  The 
particular  one  chosen  cannot  be  predicted,  since  in  the  preceding  chaotic  interval  the  solution  path 
wanders  irregularly  over  the  permitted  region,  and  selects  the  positive  or  negative  well  in  an 
apparently  random  manner. 

Publications  [11-14]  present  the  main  results  for  the  2DoF  model.  They  should  be  compared 
with  those  of  a  uniform  beam  on  the  one  hand,  and  the  1  DoF  model  on  the  other.  Paper  [11] 
includes  illustrations  of  uniform  beam  time  histories  for  comparison  with  those  of  the  2DoF  model, 
for  which  the  quasiperiodic  and  chaotic  behaviors  are  identified  by  characteristic  phase  plane 
diagrams  and  plots  of  power  spectral  density  and  Lyapunov  characteristic  exponents.  The  use  of  the 
energy  approach  is  detailed  in  publications  [  1 2]  and  [13],  with  pictures  of  strain  energy  surfaces  for 
sevei^  cases,  and  comparisons  between  the  characteristic  diagrams  (final  deflection  plotted  as 
ftmction  of  the  load  pulse  parameter)  of  the  single  and  two  degree  of  freedom  models.  The  orderly 
pattern  in  the  former  case  contrasts  with  chaotic  irregularity  in  the  latter  [12].  A  major  contribution 
of  [13]  is  a  pair  of  diagrams  showing  the  correlation  between  the  sign  of  the  largest  Lyapunov 


exponent  and  the  position  of  the  total  energy  plane  relative  to  the  eiiergies  at  the  five  stationary 
pomts  of  the  strain  energy  surface.  Also  briefly  discussed  in  [  1 3]  is  an  ^tmiative  problem 
definition,  in  which  the  four  plastic  strains  are  regarded  as  independently  chosen  parameters,  and 
the  reponse  behavior  is  studi^  for  various  initial  displacement  and  velocity  conditions.  I^liminary 
calculations  show  the  sensitivity  of  the  character  of  the  response,  which  may  be  either 
quasiperiodic  or  chaotic  at  a  fix^  energy  level,  to  the  orientation  of  the  initial  displacement  vector. 

Laboratory  experiments  were  a  fundamental  part  of  the  investigation.  A  summary  of  our 
experimental  resulte  is  given  in  [  15].  Anomalous  rinal  deflections  were  obtained  in  tests  on 
specimens  of  aluminum  alloy  6061 -T6  with  fully  fixed  ends.  These  were  somewhat  more  difficult 
to  obtain  than  anticipated.  TWs  was  probably  bemuse  effects  of  damping  were  stronger  than 
expected;  calculations  [9]  had  shown  that  anomalous  behaviors  tended  to  be  suppressed  by  energy 
losses  in  damping.  Some  progress  was  made  toward  finding  a  method  of  estimating  the 
circumstances  in  which  anonialous  behavior  may  be  expect^. 
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